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Phyllotaxis, the search for the most homogeneous and dense organizations of small disks inside 
a large circular domain, was first developed to analyze arrangements of leaves or florets in plants. 
Then it has become an object of study not only in botany, but also in mathematics, computer 
simulations and physics. Although the mathematical solution is now well known, an algorithm 
setting out the centers of the small disks on a Format spiral, the very nature of this organization 
and its properties of symmetry remain to be examined. The purpose of this paper is to describe 
a phyllotactic organization of points through its Voronoi cells and Delaunay triangulation and to 
refer to the concept of defects developed in condensed matter physics. The topological constraint 
of circular symmetry introduces an original inflation-deflation symmetry taking the place of the 
translational and rotational symmetries of classical crystallography. 
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INTRODUCTION 



The densest organization for a packing of small disks on an infinite plane is obtained when their centers are at the 
nodes of the triangular tiling of a hexagonal crystalline lattice, they all occupy the same area on this plane. If the 
disks are to be organized within a finite compact domain of the plane, this solution stays valid only if the boundaries 
of this domain are aligned along crystallographic directions. In the case of domain with a circular border, a constant 
area per center can only be obtained when the centers are regularly placed on the spiral drawn by the algorithm of 
phyllotaxis. 

Phyllotaxis in order to tile surfaces as a circular domain of the plane, a cylinder or a sphere leads to the building of 
distributions of points having the best homogeneity and isotropy possible [lJ-Q - This solution appears in plant growth 
as being the only one compatible with a sequential growth within the frame of a self organizing process submitted to 
geometrical constraints [J-Q. This has been noticed by D'Arcy Thompson^ who write: "...and not the least curious 
feature of the case is the limited, even the small number of possible arrangements which we observe and recognise." 

Also for the optimization of homogeneity and isotropy, it appears in the formation of Benard-Marangoni convection 
cells in cylindrical containers [H Q , the organization of ferrofluid droplets falling down in silicone oil in presence of an 



inhomogeneous magnetic field with cylindrical symmetry [28( and that of air bubbles on a circular water surface [11 
Finally, spherical phyllotaxis was used to estimate the Earth coverage of satellite constellations [12 1. 

A phyllotactic structure on a surface is a set of points which define position of physical objects as disks, convection 
cells, florets or others. We go into this matter describing a phyllotactic organization of points of the plane through 
its Voronoi cells and Delaunay triangulation and making call to the concept of defects developed in condensed matter 
physics. The topological constraint of circular symmetry introduces an original inflation-deflation symmetry taking 
the place of the translational and rotational symmetries of classical crystallography. 

The phyllotactic structure is represented in Figs. [2[ [31 Outside a core, it consists of grains of hexagonal cells (in 
red) that are concentric circular rings, bounded and separated by circular grain boundaries made of n-i heptagonal 
cells (inner), n\ hexagonal cells (middle) and 112 pentagonal cells (outer). The two numbers n\ and ni are successive 
Fibonacci numbers. These grain boundaries serve as natural boundaries for our optimal packing problem. Their 
properties which are presented here therefore solved the optimal packing problem topologically and metrically. 

At the origin of our work are questions raised by the structure of collagen fibrils. They can be considered as dense 
packings with circular sections of more or less parallel rods, the so-called triple helices, whose lateral organization 
was the object of several investigations through X rays scattering studies. These observations could not agree with a 
purely hexagonal lattice and revealed the presence of an important degree of disorder. A lattice built from a distorted 
hexagonal lattice with a large unit cell, a quasicrystalline structure and paracrystallinc fluctuations in a hexagonal 
lattice were then proposed. Those propositions are however difficult to conciliate with a twisted organization of the 
fibers, not compatible with the propagation of any long range lateral order, and the roundness of their normal section, 
as such orders should lead to facetted sections. As some electron micrographs also suggest a spiral assembly, we 
recently examined a model in which the sections of the triple helices would be organized according to a phyllotactic 
pattern [l3| . Even if no definitive answer can be given, main features of the X ray diffraction patterns modeled using 
phyllotaxis are compatible with experimental results. We are actually investigating this problem. 



PHYLLOTAXIS 



The generative Fermat spiral 

Sites, the points at center of physical objects are lying on a spiral (figurc[T|) define by the equations relating cartesian 
coordinates x, y with polar coordinates p(s) = a^fs and 0{s) = 2tt\s: 

x(p,8) = a^/s cos(27rAs) 

y(p, 0) = a*Js sin(27r As) (1) 

where a is a parameter defining the metric scale, A an important parameter which will be discuss below and s the 
parametrization of the curve. Sites, indexed by integers s, are placed on the curve, so that the azimuth between two 
successive sites varies by 2-7rA. The important result of the mathematical studies is that an homogeneous distribution 
of sites is only obtained with A = 1/t the inverse of the golden ratio r = (1 + \/5)/2, or A = 1/r 2 . The golden ratio 
is defined by 1/r 2 — 1/t — 1 = so that this second value gives the same result as the two corresponding angles have 
the same absolute value modulo 2n for integer s. 



FIG. 1. A Fermat spiral with 50 points on it, successively placed at angle 6 = 2nsX, with s = to 49. 



Two characteristic properties of this generative spiral are important in order to have the best covering. First, 
the added area with any new site must be a constant, for any A: this results of p(s) oc y^s, but this choice is not 
sufficient to impose a local isotropy of the distribution of sites. It is the second properties of the generative spiral 
which gives the best local isotropy of the area associated to the best uniform density: the choice of A = 1/r. With a 
rational A, sites are aligned on straight radial spokes and with irrational A different from 1/r (or of a noble number) 
they form a few long spiralling spokes, so in these two cases sites are gathered on lines jil. fl4|. Uniformity is indeed 
associated to the property of "noble number" to which t and r 2 belong, to have only 1 in the tail of their continued 
fraction expansion. Then these numbers are approximated by successive truncations of their continuous fraction which 



converge smoothly |15| . This point is discussed in [16j , which defines the best uniformity by the property of having 



shape and area of Voronoi cells as independent of s as possible, a context-free inflatable structure; see also the so 
called "shape invariance" in fl8j | . 

We presents here structures in which sites define local domains of the same area, but there are also examples where 
the area change with the radial position, like it is in a daisy where external florets are larger than those of the core. 



These organizations are obtained with other kind of spiral equations, for instance logarithmic spiral [ 1 81 ] . Nevertheless 
conformal transformations but also shearing allow to explore these different phyllotaxis. This is discussed at the end 
of this text. 



Voronoi tiling 

The geometry of a set of sites can be very well revealed using a Voronoi decomposition of the space (figures [2] 
and [3]). In 2D, each site is surrounded by a polygonal cell which is the locus of all the points of the plane closer to 
the considered site than to all other sites. In particular the Voronoi decomposition define strictly neighbors of a site 
whose cell share an edge with the cell of the considered neighboring site. Building a Voronoi decomposition for a 
finite set of sites, we face the difficulty of defining neighbors close to the border of the set. We have solved this, by 
doing the decomposition for N sites embedded in a slightly larger set. With the generic hypothesis of a coordination 



FIG. 2. Voronoi decomposition of a 2500 sites model. Voronoi cells define the number of neighbors: 6 for hexagonal cells in 
red, five for pentagonal cells in blue and seven for heptagonal cells in green. Small black points are for sites over 2500 added 
in order to avoid difficulties encountered to define Voronoi cells at the border of a finite set. 

number c = 3 for Voronoi cell vertices, topological constraints are imposed by the Euler relation F — E + V = x, with 
the Euler-Poincare characteristics x = f° r the infinite plane or for any closed surface without Gaussian curvature, 
like a torus. These constraints are such that the average number of edges of Voronoi cells to be strictly 6. But a 
phyllotactic pattern is a finite set with a circular boundary. For a compact finite part of the plane, the Euler relation 
V — E + F — x with the Euler-Poincare characteristic x = 1 has to be considered. 

In phyllotaxis there are hexagonal cells, pentagonal cells and heptagonal cells only. Thus, a tiled plane surface en- 
closed inside a circle, in the limit of an infinite radius must have the same number of pentagonal and heptagonal cells. It 
appears that the three kind of cells are organized by blocks. Outside a core, it consists of grains of hexagonal cells that 
are concentric circular rings, bounded and separated by circular grain boundaries (f u -i, fu-2, fu-i) made of f u -i hep- 
tagonal cells (inner), / u _ 2 hexagonal cells (middle) and f u -i pentagonal cells (outer). The f u are Fibonacci numbers 
defined by relation f u = + f u -% with /i = 1 and f% = 1 that is the sequence 1,1,2,3,5,8, 13, 21, 34, 55, 89, 144, .... 
These grain boundaries serve as natural boundaries for our optimal packing problem. Outwards packing begins with 
the first complete grain boundary (13, 8, 13) with 13 heptagons, 8 hexagons and 13 pentagons. The core disk is 
bounded by the 8 pentagons of the (first) incomplete grain boundary (3, 5, 8). It has 3 heptagons instead of the full 8. 
The additional pentagons (nearly) fulfil the topological requirement that a tiled disk should have a topological charge 
of 6 (i.e. 6 additional pentagons, a sphere has a topological charge 12) 0, [l9[. 

To resume, the hexagons are either in crystalline large grains (in a topological sense) or form narrow rings between 
heptagonal and pentagonal rings. The number of hexagons in the large rings are ...34, 134,422, 1221,3384, ... as we 
shall sec (tabic 1). 

RINGS AND DEFECTS 
Grain boundaries and defects 

In disordered material there are two kinds of disorder refereing to a perfect crystal. First is the metric disorder 
such that relations of neighborhood are conserved but with change for the length between neighbors, then topological 
disorder which imply defects like dislocations, disclinations breaking locally these relations. The decomposition of the 
phyllotactic pattern into rings gives good examples of these two kinds of disorder (figures [2] and [3]). The large rings of 
hexagonal cells have metric distortions but without defect. They can be seen as a ribbon cut into a perfect hexagonal 
crystal and then flattened onto a flat circular ring gluing its two ends, an operation which introduces change of lengths 
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TABLE I. Table of cell types in the successive rings or annuli. Neighbor separations Ss are Fibonacci number, all the same for 
a given ring, except at the beginning. A label u for large annuli corresponds to the Fibonacci number f u giving the medium 
separation (/„ in the set of fu, fu+i) positive separations). Rings corresponding to defects are marked with ||. 



u 


Cell type 


number of cells 


s from 


to 


neighbor separations Ss 




pentagon 


2 





1 


1,2,3,4,5 or -1,2,3,5,8 




hexagon 


1 


2 


2 


-2,2,3,5,8,13 




heptagon 


3 


3 


5 


(-3,-2,2) or (-4,-3,-2) or (-5,-3,-2),3,5,8,13 




hexagon 


1 


6 


6 


-5,-3,3,5,8,13 




pentagon 


2 


7 


8 


-5,-3,5,8,13 




hexagon 


1 


9 


9 




II 


hexagon 


5 


10 


14 


-8,-5,5,8,13,21 




heptagon 


3 


15 


17 


-13,-8,-5,5,8,13,21 




hexagon 


5 


18 


22 


-13,-8,-5,8,13,21 




pentagon 


8 


23 


30 


-13,-8,8,13,21 


7 


hexagon 


2 


31 


32 


-21,-13,-8,8,13,21 


|| 


heptagon 


13 


33 


45 


-21,-13,-8,8,13,21,34 




hexagon 


8 


46 


53 


-21,-13,-8,13,21,34 




pentagon 


13 


54 


66 


-21,-13,13,21,34 


8 


hexagon 


34 


67 


100 


-34,-21,-13,13,21,34 




heptagon 


21 


101 


121 


-34,-21,-13,13,21,34,55 




hexagon 


13 


122 


134 


-34,-21,-13,21,34,55 


|| 


pentagon 


21 


135 


155 


-34,-21,21,34,55 


9 


hexagon 


134 


156 


289 


-55,-34,-21,21,34,55 


|| 


heptagon 


34 


290 


323 


-55,-34,-21,21,34,55,89 


|| 

1 1 


hexagon 


21 


324 


344 


-55,-34,-21,34,55,89 


|| 


pentagon 


34 


345 


378 


-55,-34,34,55,89 


10 


hexagon 


422 


379 


800 


-89,-55,-34,34,55,89 




heptagon 


55 


801 


855 


-89,-55,-34,34,55,89,144 




hexagon 


34 


856 


889 


-89,-55,-34,55,89,144 


|| 


pentagon 


55 


890 


944 


-89,-55,55,89,144 


11 


hexagon 


1221 


945 


2165 


-144,-89,-55,55,89,144 


|| 

1 1 


heptagon 


89 


2166 


2254 


-144,-89,-55,55,89,144,233 


II 


hexagon 


55 


2255 


2309 


-144,-89,-55,89,144,233 




pentagon 


89 


2310 


2398 


-144,-89,89,144,233 


12 


hexagon 


3384 


2399 


5782 


-233,-144,-89,89,144,233 


II 


heptagon 


144 


5783 


5926 


-233,-144,-89,89,144,233,377 




hexagon 


89 


5927 


6015 


-233,-144,-89,144,233,377 


II 


pentagon 


144 


6016 


6159 


-233,-144,144,233,377 


13 


hexagon 


9167 


6160 


15326 


-377,-233,-144,144,233,377 



(on a cylinder, it would be possible to wrap it without distortions). With reference to a perfect hexagonal 2D crystal, 
pentagons and heptagons are disclinations of opposite weights, but a dipole of opposite disclinations is a dislocation as 
shown on figure [4j In the Voronoi decomposition pentagons and heptagons form dipoles which are dislocations in the 
structure. Several dislocations gathered along a line defined a grain boundary between disoriented grains. Hexagonal 
cells which are between two dipoles inside the dipoles ring, control the distance between dislocations, and then the 
disorientation between the two separated large hexagonal rings which are called crystalline grains even if they have 
metric disorder but without topological disorder. Dislocations are an important factor to kept the best homogeneity 
of the structure. 
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FIG. 3. Voronoi decomposition of a 6200 sites model. 




FIG. 4. A dislocation in an hexagonal crystal. The figure is drawn with regular polygons, so leading to the appearance of a 
stair. In 2D the structure is flatten by projection on a plane. 

Parastichies 

Parastichy is a botanical term used to design spiral lines observed on some plants such as composed flowers, cactuses 
and pine-cones. Looking at figures [2] and |3] such lines appear, even if they are not easy to follow all along. Some 
segments of parastichies appear on figures [5] and |6] The purpose is now to describe them precisely. Points s and s + 1 
placed on the generative spiral are separated by long distances (at least for not too small s). For visible parastichies 
we are concerned by points which are neighbors, as defined by Voronoi cells, and search for the separation Ss between 
two such points labeled s and s + Ss. Two points with polar coordinates (p, 9) are close if their p and their 9 are 
close. So ay/s must be close to a\J s + Ss, implying a small 8s /(2s). The two azimuthal separations must follow 
2tt\s — 2ix\(s + Ss) ~ modulo 2n. It is impossible to have in this equation a strict equality to zero because A, the 
inverse of the golden ratio, is an irrational number but it can be approximated by a rational number A' = f u -i/fu, 
the ratio of two consecutive Fibonacci number. The equation for the azimuthal separation, with this approximation 
reduces to 27rA's — 2irX'(s + Ss) = which is solved exactly by 27r/„_i = modulo 2-7T for 8s — f u . The approximation 
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FIG. 5. Family of lines formed with parastichy segments, all rotating with the same orientation, corresponding to the smallest 
separation Ss given by a Fibonacci number /„ with an even it (for conservation of the rotation orientation). The color code 
is the same as on figure 2 and 3: red for 6 neighbors, blue for 5, green for 7. A new line originates when a new Fibonacci 
number, with these properties (the smallest with even u), appears in the list of separations with neighbors. Except for very 
small s, this happen on sites with pentagonal cells. These sites are gathered in one over two rings of pentagonal cells, such that 
in heptagon-pentagon dipole the new Fibonacci number is not in the list of separation for the site with the heptagonal cell, 
but the line is oriented as the dipole. All these lines play the role of reticular lines found in 2D crystals and clearly show the 
dislocations as an insertion of extra line. 

for A is better for large u, but as Ss /s must be small, first neighbors correspond to the best compromise between 
radial and azimuthal contribution. 

Inside a ring of Voronoi cells of the same type, all neighbors are separated using the same set of Fibonacci numbers 
(except for s < 7). For instance for the ring of hexagonal cells between s — 67 and s = 100 all six neighbors are at 
s + Ss with Ss in the set {-34, -21, -13, 13, 21, 34} of Fibonacci number f g = 34, / 8 = 21, f 7 = 13 (table 1). 

Passing through a site of large hexagonal rings there are three parastichies defined by three consecutive Fibonacci 
numbers f u -i, /«, f u +i corresponding to three positive Ss separations. It is convenient to label large hexagonal rings 
by u such that f u correspond to the medium separation in this ring. So in the example with 67 < s < 100 the label 
is u = 8 corresponding to the Fibonacci number /g = 21 in the set (13, 21, 34)(see table 1). 

Equation of parastichies 

In a grain three visible parastichies are defined by three Fibonacci numbers f u —i, fu, fu+i, but now we extend the 
notion of parastichies to a whole family of curves joining sites separated by any Fibonacci number. So a parastichy 
is defined by a given Fibonacci number f u . We are showing that the parastichies equation has the form: 

x(p,6) = a^cos(0 n + {-l) u+1 2ir\ u s/f u ) 

y(p,9) = ayGsm(9 n + (~l) u+1 2ir\ u s/ f u ) (2) 

The phase shift 9 n is adjusted in order to have the curve going through the site s, so n is numbering all parastichies of a 
same family. Here we consider s as a continuous variable, in order to define a curve, but we assure that the curve pass 
through sites for integer s values. Consider the two points s and s + f u lying on a parastichies and on the generative 
spiral. On these curves their azimuthal separation is 2n\f u . This number which is largely greater than 2tt can be 
reduced, modulo 27r, to a small number between —n and 7r given by 27r(/„A — f u -i) because A ~ A' = f u —i/f u - The 
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FIG. 6. Family of lines formed with parastichy segments, corresponding to the largest separation 5s given by a Fibonacci 
number f u . Except for very small s, new lines originate on sites with pentagonal cells. The number of cell on the border is the 
number of these parastichies cut by it. 

ratio A" appears due to a property of power of the golden ratio: r" = f u T+f u -i, if A is the inverse of the golden ratio it 
followsQA" = (-l) u+1 (/ u A-/ u _i). Introduc ing the expression for A", this azimuthal separation for s varying from s to 
s + fu is (-l)" +1 27rA u . Consequently the Fermat spiral {aVicos((-l)"+ 1 27rA tl s// u + 6»„), ayisin((-l) tl+1 27rA u s//„ + 
9 n )} which rotates slowly compare to the generative spiral goes through the two points (with a good choice of the 
constant 6 n ) when s vary continuously. Notice that the scaling factor a is the same for all spirals including the 
generative spiral. A factor (— changes the orientation of the spirals from even u to odd u. This can also be put 
in relation with the fact that the approximation A ~ fu-i/ fu is an over estimation for odd u and an under estimation 
for even u. 

A generalization of parastichies: reticular lines of the phyllotaxis 

In botany, parastichies are visible spirals of the structure and so are those joining close neighbor sites. In hexagonal 
rings there are three such lines characterized by three consecutive Fibonacci number, they arc clearly visible when 
their characteristic Fibonacci number appears in table 1. Nevertheless for large s even if they are less visible they 
always exist as continuous spirals joining sites not necessarily first neighbors, but second, third... neighbors. So to 
each possible Fibonacci numbers /„ corresponds a family of f u identical spirals joining sites with a separation Ss = f u . 
There are f u in the family, because if one goes through the site s and also through the site s + f Ul there is one through 
the site s + 1 which is different, and one through s + 2 and so on up to s + f u — 1. 

The first example is the family corresponding to f\ = 1 with only one member, the generative spiral joining site s 
and s + 1. The next example is given by f% — 1 also joining sites s and s + 1. These are the two first families with a 
single spiral in each. In fact these two different spirals correspond to the two possible choices of A given by 1/r and 
1 /t 2 giving the same structure with two spirals with reversed rotation orientations 0. Exactly like for reticular lines 
in a 2D crystals all sites arc on a member of each family and lines of the same family arc in some sense parallel lines: 
in crystals they are equivalent by translation, here they are equivalent by rotation. In 2D crystals, a reticular family 



1 These two relations can be shown easily by induction considering that they are true then calculating t" +1 or A u + J using t 2 = 1 + r or 
A 2 = 1 — A. 

2 It is very usual in literature to use A = 1/t 2 to describe phyllotaxis. This choice corresponds to the smaller divergence angle. The other 
choice, done here is formally simpler and strongly supported by this property: spirals are defined by the rank u of Fibonacci number; 
their speeds of rotation around the origin decrease with u and we choose for generative spiral the first one given by u = 1. 
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FIG. 7. Variation of the area of Voronoi cells for points between s = 34 and 2166. Colors correspond to the type of Voronoi 
cells. The area are close to the mean value 7r with a = 1, strong variations appear in rings of defects. Nevertheless these strong 
variations decrease as the area tends toward it for large s. 



area 




is characterized by two Miller indices but the family of spiral is just characterized by a Fibonacci number. To resume, 
the number of spirals of the same family f u is the Fibonacci number /„. As the golden ratio r is approximated by 
the ratio of successive Fibonacci number converging toward it alternatively up and down, this result in an alternative 
rotation orientation for the spirals depending on the parity of u in /„. 



Euler relation 



As already said above, the Euler relation V — E + F = \ imposes the average number of edges per Voronoi cell 
to be 6 on the infinite plane for which x — 0. However, for a compact finite part of the plane, as is the circular 
domain of a phyllotactic pattern, x — 1 so that the number of cells with five or seven edges are not equal. We 
introduce the number of faces F — F 5 + F e + F 7 , F p being the number of faces with p edges, and the number of edges 
E = (5F 5 + 6F e + 7F 7 )/2 + (m 2 + m 3 )/2 or in order to eliminate V, E = 3(V - m 2 - m 3 )/2 + (2m 2 + 3m 3 )/2. The 
number m2 and m 3 are the number of vertices of the cells on the border with coordinations c r = 2 or 3. This gives 
(ni2 — m 3 ) + (i*5 — Fj) = 6. Two particular solutions of this equation are interesting. First, following what happen 
on the infinite plane, suppose that F$ = F? so there is mi = 6 + m 3 , an excess of six sites with coordination 2 on the 
border. A simple example is obtained with a finite piece of an hexagonal tiling cut along edges of hexagonal cells ( so 
with F5 = Fj — 0). Then this structure has an irregular boundary depending on the choice of hexagons kept inside it 
respecting 7712 = 6 + to 3 for the ui2 sites of the boundary belonging to only one hexagon, and the m 3 sites belonging 
to two hexagons. The other particular solution, corresponds to the phyllotaxis: the boundary is close to a circle with 
TO2 = m 3 , an alternative cquipartition of coordination 2 and 3 on the border. Notice that in a perfect hexagonal tiling 
cut along edges in order to divide it in two pieces with sites on the cut alternatively on one or two hexagons (c r = 2 
or 3), the cutting line oscillates around a straight line, so we can consider such border as some kind of "topological 
geodesic" 0. In this case, solving the equation, gives F 5 — F 7 = 6 corresponding to an excess of six pentagons in 
the core of the structure. Evidently a circular border with these properties (mi = m 3 forcing F5 — F 7 = 6) must be 
inside large hexagonal domain and not inside rings of defect: otherwise there would have an excess of heptagon. A 
consequence of m 2 = m 3 on the border, is also that the number of cells on the border (having outside edges) does not 
varies continuously with the total number of site in the phyllotactic pattern: this number remains constant in large 
hexagonal rings, but otherwise changes rapidly when the boundary enter inside defect rings. This number is equal 
to the number of parastichies corresponding to the highest separation between neighbors, so if the border is inside 
non-defective domain it is equal to Fibonacci numbers. 
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FIG. 8. Delaunay triangulation in a small region around the point s = 15326 which is the last point in a wide six neighbors 
ring. The center of the phyllotactic disk is far away on the left. Red points are common to 6 triangles, blue points to 7 triangles 
and green points to 5 triangles. Often it seems that two triangles form a square, but in fact in this sheared square the diagonal 
which divide it into triangles is slightly shorter than the other one. There are two domains differentiated by the diagonals 
orientation. In these domains all points are common to 6 triangles. Between these two domains the defect ring appears as a 
place where some diagonal orientations are flipped, leading to 7,5 or eventually 6 triangles around a point. Going away from 
the defect domain, triangles become more and more regular. 

Rings and defects: a relationship with crystals 

Exactly like for reticular lines in a 2D crystals all sites are on a member of each family and lines of the same family 
are in some sense parallel lines: in crystals they are equivalent by translation, here they are equivalent by rotation. 
In 2D crystals, a reticular family is characterized by two Miller indices but the family of spiral is just characterized 
by a Fibonacci number. To resume, the number of spirals of the same family f u is the Fibonacci number /„. As the 
golden ratio r is approximated by the ratio of successive Fibonacci number converging toward it alternatively up and 
down, this result in an alternative rotation orientation for the spirals depending on the parity of u in f u . 

There are relationships between crystals with defects and phyllotactic pattern, but in a next section relations to 
quasicrystals are given. 

METRIC PROPERTIES 
Packing density and fluctuations 

As the circle of radius p and area irp 2 = ira 2 s contains s points, the area added per point has the value 7ra 2 .The 
area of Voronoi cells oscillates close to this value for small s and then converges towards it (fig. [7]). 

Influence of site distances on shapes of Voronoi cells 

In the Voronoi decomposition of a phyllotactic pattern there are only hexagons, pentagons and heptagons, nev- 
ertheless, looking rapidly at figure 3 some squares seem to appear, mainly in rings related to dislocations. In fact 
these are not squares, but Voronoi cells which could be described as a square with one corner slightly truncated (for 
pentagon), or with two truncated corners (for hexagon), or three truncated corners (for heptagon). This is to relate 
to the fact that the two parastichies corresponding to the two first Fibonacci separations are close to be orthogonal 
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FIG. 9. First neighbor distances between sites s and s + Ss. Blue, for the interval Ss equal to the smaller positive Fibonacci 
number in the list of table 1. Green, for the next interval, red for the third positive interval and purple for the last one (occurring 
only if the Voronoi cell is an heptagon). So green, blue and red correspond to distance along the three visible parastichies. Each 
continuous curve corresponds to a given Fibonacci number that appears in different annuli. For instance, /n = 89 appears 
between s = 290 and 5926 leading to a continuous curve which is successively purple, red, green and blue. Each lower crossing 
of two curves corresponds to a grain boundary, upper crossing is in the middle of a hexagonal annulus. The scaling parameter 
in equation (1) is a = 1, so the distances are approximatively in the range [1.676, 2.506] with a mean value 1.903. 
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in rings of defects. Oppositely in large rings of hexagonal cells, the core of the ring contains less irregular hexagonal 
cells. 

Figure [5] shows how Delaunay triangles are organized near the rings of defects In these regions corresponding 
Voronoi cells are close to be square, but with corner slightly truncated. Effectively Voronoi cells vertices are centers 
of circumcised circles of Delaunay triangles, so if two such triangles look like a square divided by a diagonal the two 
centers are close to the diagonal leading to a small edge of the Voronoi cell. If two corners are truncated the cell is an 
hexagon. In large hexagonal rings, the Voronoi cells are rhombus with two opposite truncated corners, so the Voronoi 
cells are close to regular hexagons. In the narrow hexagonal rings, inside defect rings, the cells look like squares with 
two close truncated corners, so are very distorted hexagons. 

Even if Delaunay and Voronoi decomposition give exactly the same information some properties are more easily 
observed with one or the other method. For instance angles between parastichies are very clear on the figure [8] because 
edges of triangles are along parastichies. It is clear in this region of defects that two families of parastichies are close 
to be orthogonal with a third family making a 7r/4 angle with the others. 



Distances between first neighbors 

The behavior of distances between first neighbor sites is shown on figure |9] In domains where Voronoi cells are 
slightly deformed squares, their area is approximatively ira 2 so that the distance between points, along a square edge, 
is s/tt ~ 1.772 (with a = 1) and along a diagonal, is s/2tt ~ 2.506. When Voronoi cells are more clearly hexagons two 
distances are close to J ^ — 2.053. The behavior of the distance between the two points defined by s and s + f u as 
a function of s is obtained using relation given by Yeatts j20|: 

d u {s) = fu{— + J} ) 1 , (3) 



Delaunay triangulation is such that sites are associated to cover the surface by triangles. It is related to Voronoi decomposition, because 
vertices of Voronoi cells are centers of the triangles, Delaunay and Voronoi decompositions arc dual. 
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this relation perfectly fits the numerical values given on the figure |4j It can be checked that minimal values are 



all close to and converging toward \J~^ — 1-67. 

All distances are in the range [y'-^L, \/2tt}, whatever will be s and u. This property, specific to the choice of A = 1/r 
is very important to ensure the best packing homogeneity 0. 



GRAIN BOUNDARIES 



Table 1 shows lower and upper bounds of domains containing the same kind of Voronoi cells. These values are 
obtained from numerical observations of a phyllotaxis. Nevertheless it would be interesting to have an analytic relation 
giving these bounds. Clearly having the lower bound for hcptagonal or pentagonal rings the upper bound is obtained 
adding a Fibonacci number minus one. This is also true for narrow hexagonal domain, in the same region, using the 
previous Fibonacci number. But there is no clear relation for wide hexagonal rings. We introduce a semi-empiric 
relation which has been tested up to s = 280783. 

As in a ring of defects Voronoi cells are close to be squares, we suppose that sites are with a good approximation 
on a square lattice with parameter d = yfn. We set a = 1 in the following. Figure ITUl a shows a representation of 
the defect ring at 101 < s < 155 in which Voronoi cells are 21 heptagons, 13 hexagons and 21 pentagons. It appears 
straight in this representation, but in the phyllotaxis it is refolded into a ring. The defect ring appears as a strip of 
length given by the vector (21, 34) in the square lattice and a width given by (1,-1), so that the number of points in 
the strip is 55, given by the determinant constructed with the two vectors. This strip is separated into three thinner 
strips of same length, but with widths defined by vector (0, —1), (1, 1) and (0, —1) respectively. The first one contains 
the 21 points with heptagonal cells, the second contains 13 points with hexagonal cells and the third contains 21 
points with pentagonal cells 0. In the phyllotactic pattern, the 55-points strip is folded into a flat ring that bounds a 
disk, by identification of the two small sides at the price of metric distortions. We estimate the bounds of the strip 
using two properties: the length of the strip and the number of points in the disk that is enclosed with radius s/s. The 
length is the module of the vector (21, 34) which is L = V21 2 + 34 2 y / 7r using d = y/% as length unit. Folding the strip 
into a ring, this length is the perimeter of a circle. We can estimate the number of points of the phyllotactic pattern 
in the disk enclosed by this circle using the definition of the generative spiral. This number is Sf, = 21 ^ 4 ~ 127.085. 

Now we compare this to the known value for upper and lower bounds in this rings of defects. When the strip is 
folded into a ring, the internal side of the ring is compressed whereas the external side is expended, so Sf, is collimated 
by the narrow ring of hexagons. The evaluated value Sb falls between these bounds (122, 134). 

From this estimation we deduce a formula which gives bounds of all hexagonal rings appearing inside defects regions: 

/ i r ^ j~ u 1 /u+/u+l n ,r/u— 1 + 1 . /a + /n+ln 

( round [ _ + _ ], r0 und[ g + An ]) 

It is possible to write this relation in a more compact way using the identity / 2 + / 2 +1 = fa u +i , which can be shown 
by inductive reasoning on v using f u + v = f v +ifu + fvfu-i- Then we have: 

( r o U n d [ g + ^T^ TOmid ^ ~ 2 + ( 5 ) 

Here f u is the Fibonacci number of order u and the function round take the integer part of an irrational number. The 
example of the figure [TOl corresponds to /§ = 21. In this relation the two bounds are such that the number of point 
in the ring is / u _i as wanted. This formula has been checked to give exact value (except a shift of one for u = 11) 
up to the bounds (280174, 280783) obtained with u = 16. Using the fact that rings involved in defect domains have 
number of points given by Fibonacci number it is easy to get all other bounds as given in table 1. 



4 These limits remain the same in phyllotactic tilings on the sphere and on the hyperbolic plane, providing the same it area for each site. 
This will be given in a further article on phyllotaxis in non-Euclidean geometry. 

5 Notice that this strip is very similar to the construction of approximants of quasicrystal using the cut and projection method. It is one 
of the reasons which relate phyllotaxis and quasicrystal. 
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FIG. 10. Strip cut in a square lattice (a). This strip is a good approximation of the ring of defects between points s = 101 and 
s — 155 containing 55 points. Color of the points correspond to the type of their Voronoi cells. Point s = 101 is the first green 
point down, then points are numbered, increasing by 21 going up or decreasing by 34 going right. This strip can be divided into 
three strips of heptagons, hexagons and pentagons, (b and c) Continuation of the grain toward the other boundary on the same 
square grid {(x, y) : x — 1 ■<-> s + 34, y — 1 ■<-> s + 21} which can be considered as folded on a cylinder with axis perpendicular 
to the strip. This grain (in grey) is limited by an inner border (b) and an outer border (c). Blue points (pentagons) in (b) are 
those of (a), orange points belong to the large grain. The grain is bordered in (b) by a Fibonacci chain defined by the lower 
part of the grey domain. There are 21 stairs in it containing 21 blue points and 34 orange points (hexagons of the large grain). 
The other border of the grain (in (c)) is geometrically identical, but it contains 34 green points (heptagons of the next ring of 
defects) and 21 red points (hexagons of this ring of defects). This next ring of defects from s = 290 to s = 378 contains 89 
points. The beginning of this ring is expanded in (d), with two Fibonacci chains, one which close the grain from s = 156 to 
s = 289 and the other which begin the next grain, from s = 379 to s — 800. The two Fibonacci chains could be coded with long 
(purple) and short (blue) segments forming a "stair case". Stairs of the first chain appear with n/2 angles, but for the other 
angles are 7r/4. In fact the strip should be deflated (stretched along by a rational approximant of r and compressed across by 
an approximant of 1/r, that is a Poisson shear). Introducing progressively the Poisson shear, instead of a part of a cylinder 
the substrate of the grain and its two boundaries is a part of a surface with a negative constant Gaussian curvature around the 
same axis. This Poisson shear change the rhombic cell darken in (d) into a square. Refer also to figures 11 and 12. 




There is an other way to look at this problem. It is to search s for points where two parastichies are very close 
to be orthogonal. This relation which is deduced from results given by Yeatts [2(| is: s = abs[ ] where 7„, 

related to approximants of r, is 7„ = 2ir(\f u — f u -i)- This relation give real values for s very close to ■ /2 " +1 but 
clearly different (for instance s' b ~ 127.067 compared to ~ 127.085). 

It will be a challenge to found an exact relation giving integer values for the bounds without ambiguity on the use 
of function "round" , "floor" or "ceiling" extracting integers from non integer numbers. 

The ratio between radius of two successive rings of defects, are proportional to \J J2U+1//2U-1 very close to the 
golden ratio. 
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CONFORMAL TRANSFORMATIONS AND SHEARING IN PHYLLOTAXIS 

In this section, we show how the successive grains are sheared relatively to a perfect crystalline structure, so that 
each parastichy is a stack of deformed hexagons, starting from square at one boundary, through more symmetrical 
hexagon at mid-grain, to square at the next boundary. See figure [9l where the lowest crossings correspond to square 
cells at the grain boundary with d = ^/tt (the third distance is the diagonal v2tt ~ 2.507, the maximal distance). The 
upper crossings correspond to hexagons with 27r/3 angles (but not perfectly regular) at mid-grain, with a minimal 
distance given by |QJ. 

Phyllotaxis (1) is an example of spiral lattice [2l| . the image D in the complex plane w = pcxpiO of a regular 
lattice P in a domain of the complex plane z with cartesian coordinates. In elasticity, D is the actual, deformed state 
of the material, and P is called the natural state. The mapping is through the function w(z). 

One particular mapping w(z) = w ea;p(b*.z) is conformal, ic. it is analytic, invertible, singularity-free (excluding 
a small domain around the origin of w), and thus without any defect (dislocation, disclination or grain boundary). 
It conserves the angles and is also unsheared (a consequence of the Cauchy-Riemann relations imposing analyticity) 
21 1 . 

We consider a square lattice and we define its nodes positions by complex numbers in its plane. As an example we 
consider a strip in this lattice, like that shown on the figure [10] with square vertices selected in. We have to map it on 
a ring of dipole. On this example the strip is define by a long multiple cell of the square lattice characterized by the 
two vectors = (21, 34) and bo = (1, —1). The conformal mapping transform the imaginary axis into a circle, so we 
have to choose this axis along the b 2 vector. As complex numbers are defined using orthogonal real and imaginary 
axes, we introduce the real axis along the vector bi = (34, —21). To express the argument of the mapping function we 
introduce reciprocal vectors of h± and b 2 called h\ and b£ and defined by b;.bj = Sij. Then the mapping function 
is w(z) = w exp[27r(b^ + ibjXxbi + yh 2 )] where x and y are coordinates of a point z expressed in the multiple cell. 
The mapping function reduces simply to w(z) = w exp[2Tr(x + iy)] with the scaling parameter w 

The square lattice on P includes all points with s > 101. Note the particular families of reticular lines, (i) (i) (1, 0), 
(0, —1) and (1, —1), that are mapped into parastichies, (ii) (21, 34), parallel to the grain boundary, that are mapped 
into concentric circles, and (iii) (34, —21) that are mapped into radial spokes. The conformal mapping yields a spiral 
lattice [2l[, with all reticular lines mapped into equiangular (or logarithmic) spirals, save (21,34) (concentric circles) 
and (34, —21) (radial spokes). All angles are conserved in the mapping, notably the two parastichies, images of (1, 0) 
and (0, —1) remain at right angles and keep the same angle with the circles parallel to the grain boundary and with 
the perpendicular radial spokes. Thus, figure [TU] is a conformal (literally) representation of the grain boundary, with 
square-shaped Voronoi cells. However, the density of points \dw / dz\~ 2 — |&ui|~ 2 is not uniform and the conformal 
mapping must be sheared, from the grain boundary outwards, to achieve the required uniformity. Accordingly, the 
hitherto defect-free phyllotaxis will exhibit other grain boundaries further out. Remarks: a) the mapping is enacted 
outwards, as occurs naturally in the growth of phyllotactic structures, younger florets pushing out the older ones and 
squeezing in-between at grain boundaries; b) we can say nothing for small s. In particular, on the central core of the 
phyllotaxis s < 32 that is probably the smoothest way of accommodating the core of a weak, non-singular topological 
defect that is here extra- matter [19j. 

Figure 1 101 - a is a representation of the lattice P. It consists of adjacent, non-overlapping parallel strips of equal 
thickness generated by the translation bx- Each strip in bounded by two reticular lines parallel to (21, 34) and contains 
55 points, the diagonal of the unit square goes as the thickness of the strip. The two points are separated by 8s = 55. 
The first strip is the grain boundary itself. Successive strips are mapped conformally by w(z) = woexp[27r(a; + iy)] 
into concentric rings containing 55 points each, but with thickness increasing exponentially p(s + 55) = p(s) cxp(55), 
where p — \w\. In order to have uniform density of sites in the presented phyllotaxis with p(s) — a^/s, the successive 
strips in P must be squeezed across and extended through a Poisson shear strain (figure ITOl-d) . 

RELATION WITH QUASI PERIODIC STRUCTURES 

The way dislocation dipoles are organized along circles is strongly related to ID quasicrystals. A Fibonacci ID 
quasicrystal can be obtained using an inflation-deflation rule iteratively applied to a sequence of long and short 
segments. This inflation-deflation rule is L — > L + S and S — > L. Starting simply from a short segment this specific 
rule gives a quasicrystal after an infinite number of iterations, or with a given number of iterations, a finite structure 
with a number of short and long segments given by two successive Fibonacci numbers. Consider now dipoles along 
circles in the phyllotaxis. There are isolated dipoles (an heptagon and a pentagon in contact) and pair of close dipoles. 
Consider now the rule: change isolated dipole into a pair of dipoles and change pair of dipoles into a pair and an 
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isolated dipolc. It is the rule we have going to one ring of dipoles to the next ring. So there is an inflation-deflation 
symmetry associated with radial change relating defects in this structure. This can be checked counting the number 
of pair of dipoles or of isolated dipoles which are successive Fibonacci numbers, on circles of defects. The figure [TU1 
which is similar to figures used in presentation of one dimension quasicrystal defined by cut and projection is related 
to this property. 



Successive grain boundaries 



The first complete grain boundary (13, 8, 13) marks the beginning of a regular, inflatable phyllotaxis, where the 
circular grain boundaries are words, sequences of letters L and S, where L labels a dislocation (oriented dipole 7/5) 
and S an isolated hexagon in the grain boundary. The sequences are obtained by inflation rules S — > L, L — > LS. 
Consider three successive grain boundaries, with sequences u>i_2, Wi-i and Wi. Inflation implies that the sequence 
Wi is given by concatenation of the previous two, Wi = Wi-i • Wi-2, thus w; = LSLLSLSL is a concatenation of 
Wi-i = LSLLS with Wi-2 = LSL. The word Wi has \wi\ = /, letters, d and /j_2 s. It starts with LS and 
ends with LS for i > 3 with even fi, with SL for i > 3 with odd fi. A circular word is also a palindrome. Let w T 
be the sequence w read backwards (transposed in matrix notation) . A palindrome p = p T is sequence of letters that 
reads the same either way, e.g LSLSL, hannah or eve. The word Wi can also be written as the concatenation of two 
palindromes Wi ~ pa *Pb, with w^-i = paSL,pb — SLwi-2 for i even, and = paLS,pb = LSwi-2 for i odd. 
Thus \pa\ = fi—i — 2, \pb\ = fi-2 + 2. For inflation to be consistent from one grain boundary to the next, one requires 
\pa\ > \pb\, so that the shortest complete grain boundary is (13,8, 13), a palindrome of 21 letters L and S, written 
on a circle and readable in either sense, since w = pa • Pb = Pb • PA (by circular invariance) equals w T = p^ • p T A . 

Thus, p\ = P^LSp^ 2 = p^ 2 [SL]p^[LSy A 2 for i even, and p\ = P^SLp 1 ^ = p\ 2 [LS\p l ^[SL\p l X 2 for i 
odd, with p\ = 1 



Boundary of large hexagonal rings 

A large hexagonal ring which is lying on the plane is nevertheless topologically equivalent to a perfect hexagonal 
tiling wrapped around a finite cylinder limited by two parallel circles. With such distortion of the ring, hexagonal 
cells are perfect regular hexagons and so this geometry can be seen as a multiple crystallographic cell of the hexagonal 
lattice with two identified opposite sides. The vectors which are the basis of the multiple cell entirely characterize the 
topology of its related large hexagonal ring. In order to precisely define the boundaries of a ring we will consider in 
it all pentagonal sites which are just before and all heptagonal sites just after associated with the following narrow 
hexagonal ring. With this definition of large ring boundaries the set of all large rings entirely cover continuously the 
whole structure. They are perfect topological hexagonal structures each defined by a crystallographic multiple cell 
which multiplicity is given by the number of points in it. The departure from the hexagonal coordination only results 
from the junction between successive rings. 

The inner border of a ring is a zigzag line containing pentagonal sites whose number is given by a Fibonacci number 
/u_i and hexagonal sites whose number is /„. The figure [11] gives two examples of hexagonal rings: one containing 
134 hexagons, characterized by the medium neighbor separation 5s = fg = 34 the other with 422 hexagons and 
8s = fio = 55. Consider such a ring as a hexagonal lattice, with points on vertices of regular triangles, wrapped 
around a cylinder and with the primitive cell define by two vectors it, b defining elementary translations from a 
point to its first neighbors defined by f u -i and f u . Then the ring can be seen as a multiple cell defined by two vectors 
expressed in the primitive cell by coordinates (/ u -i, — /«) and (x, y) where x and y depend on the ring width. The 

fu-l ■ 



area of the multiple cell, using the area of the primitive cell as unit is given by the determinant: 



-fu y 



It is 



the number of points in the multiple cell. 

Wrapping of the multiple cell is obtained by identification of the translation (/ u _ i, —f u ) to the identity which result 
of gluing of two sides of this cell. The two other sides give the two boundaries of the ring. This type of description 
shows that the two boundaries are equivalent with the translation (x,y). Counting the number of points in the ring 
comes to count the number of points in the cell but in the cell the two boundaries are equivalent, so both count for 
one. Then we have to add to the cell area the number of points on a border which is f u -i + fu- Following the number 
of points in the ring including pentagonal sites inside and heptagonal and hexagonal sites outside, as on the figure [TT1 



16 



First red First green 
#324 \ , #290 

1 • .-• First blue */ . * , # 



.-• First blue 

< V * > * 

F *» 
v* • > * 



« t • / 
.<■* • i * 



t 




a} 



V* . ••■.•••.* • J - 



FIG. 11. Two example of hexagonal rings: a) with 134 hexagonal sites (red), completed with 21 pentagonal sites (blue) on 
the inner border and with an outside border containing 34 heptagonal sites (green) and 21 hexagonal sites (red); b) with 422 
hexagonal sites (red), completed with 34 pentagonal sites (blue) on the inner border and with an outer border containing 55 
heptagonal sites (green) and 34 hexagonal sites (red). Inner border have blue points inside, with steps of zig-zag between 
bluepoints containing one or two red points. This is labeled I or s in a). Outer border have red points on the outside with one 
or two green points between two red points. The inner and outer sequence of I and s are the same, starting at the first "blue 
point" or the first "green point". Three families of parastichies are clearly visible, they are characterized by three Fibonacci 
numbers /«, fu+l) indicating the increase of indices jumping from a points to a neighboring point following a parastichies. 

These number are indicated on the figure (they correspond to it = 9 for a) and u = 10 for b) ) . Notice that if we orient border 
following parastichies defined by the smaller increase given by f u -i, the border in a) and b) are in reversed orientation: this is 
associated with the parity of u. 



is given by the determinant 



fu-i x + l 

-fu y + i 



In the example of figure ITD-a this is 



21 

-34 10 



giving 210 sites in the ring. For figure [TTT-b it is 



34 13 

-55 -5 



with x = — 1 and y = 9 



with x = 12 and y = —6 giving 545 site in the ring. 



from boundary to boundary 

It is convenient to refer rings by the number u in f u corresponding to their most visible parastichy. Their boundaries 
have f u + f u -i = fu+i points. How is related the outside border for u with the inside border for u + 1. The 
figure [T2l show this relation with the example u = 9 so with /„ = 34 and f u +i = 55. A chain can be described as 
a chain of long segments (two heptagonal or pentagonal sites separated by two hexagonal sites) or short segment 
(two heptagonal or pentagonal sites separated by one hexagonal site). On this example the first chain is coded 
LLSLLSLSLLSLLSLSLLSLS and the other is coded LLSLLSLSLLSLLSLSLLSLSLLSLLSLSLLSLS, but it 
is very important to notice that the chains are oriented in order to have an increase of site index, when jumping from 
their first site (first heptagon or pentagon) to the close one corresponding to the positive Fibonacci number f u -i for 
the first chain or f u for the other. Then the orientation of the two chains are reversed. Nevertheless as chains are 
palindromes using an other origin with a circular permutation, leads to symmetric representation. 

The inflation law relating the two chains is L — > LS and S — > L, as usual for Fibonacci chains but this breaks the 
symmetry It is then possible to have a law preserving the palindromic symmetry, but relating a chain defined by 
u — 1 to a chain defined by u + 2. This law which is L — > LSLSL and S — > LSL relate chains of lengthes f u and f u +3- 
This law works perfectly well if lengthes are even numbers. If lengthes are odd numbers the two first letters have to 
be put at the sequence end to respect symmetry. 
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FIG. 12. The outside border of the ring with 210 sites and the inside border of the ring with 545 sites. These boundaries are 
two chains of sites, one with heptagonal sites (green) separated by one or two hexagonal sites (red)have 55 sites, the other with 
pentagonal sites (blue) separated by hexagonal sites (red) have 89 sites. Starting from the first green site (cross) and rotating 
as indicated by arrows we can describe chains as a succession of long (I) or short (s) segments. 



It is possible to conclude that the boundaries are related by symmetry operations. Inside a given ring the inner 
and the outer boundaries are related by some kind of translation, but which appears in the plane as a spiraling 
transformation coupling a rotation and a radial expansion. Then going to the next ring there is an inflation symmetry 
associated to a reverse of orientation. 



NOBLE PHYLLOTAXIS VS. FALSE PHYLLOTAXIS 



Other values of the parameter A and their approximants 



All the properties we have discussed above are a consequence of the continued fraction expansion of A which gives, 
by truncation, convergent rational approximants When A is a noble number, {q\, q2, . . . , 1, 1, 1, . . .} in continued 
fraction expansion it is possible to extrapolate from A = 1/t = {1,1,1,1,...}. False phyllotaxis is obtained with other 
rational or irrational A, but never appear in nature. It is interesting to follow distances between neighbors (similarly 
to figure |9] given for true phyllotaxis) as a function of s and how they evolve from grain to grain. Consider the three 
distances between a site at s and neighboring site at s + Ssi as three functions of s with constant Ssi characteristic 
of each grain. These three functions have minimum values which are different, so leading to a dispersion of neighbor 
distances. Going from grain to grain so with new sets of Ssi two minimum are increasing but there is a smallest 
distance with a minimum decreasing with ,s. Voronoi cells, even if they remain hexagons in grains are more and more 
elongated with increasing s leading to a spider web aspect. This is the signature of the non- uniformity associated to 
non noble numbers. 

In the case of a true phyllotaxis we have seen in paragraph 3.2 that neighbors are gathered on f u Fermat spirals 
(parastichies) when they are separated by Ss = f u . These is because convergent approximants of A = 1/r obtained by 
truncation of it continued fraction expansion have the form / u _i/ f u . In the case of non noble numbers, the truncation 
of the continued fraction expansion of A at q u leads to principal rational approximant {A' = r u /p u ;r u ,p u € Z}. Then 
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some neighbors are separated by 5s = p u and gathered on p u similar spirals. Successive p u and r u arc related by: 

Pu = QuPu-1 +Pu-2- (6) 

In a grain consider a Delaunay triangle of close neighbors. The three edges of the triangle defined three directions of 
parastichies in the grain. The homogeneity of the grain require that these three families of parastichies, characterized 
by three index separations {5si, #S2, 8S3} are all the same in the grain. If the three separations are given in increasing 
order, the closure of edges of the triangle results in the triangular relation Ss^ = 5s\ + Ss2- If in a grain, a site s 
is neighboring a site s + p u the triangular relation and equation [6] are incompatible except if q u = 1 (as it is for all 
not to small u in true phyllotaxis) . To solve this incompatibility new intermediate approximants appear: to r u /p u 
are associated q u approximants r« /p$ with the sub-index (i) <G [1, «/«]. The integer p^} and r u ^ are defined by 
Pu+i = Pu-i +Pu and p^ +1 = Pu+i +Pu, so we can write = p u -\ and p£+\ = p u +i- Using denominators of 
these intermediate approximants it is then possible to satisfy the triangular relation with the separations. In this case 
at least one principal approximant denominator p u appear in this list, and it appears in q u + 2 successive grains. 

Notice that distances between sites separated by a Ss related to an approximant ri* /pi are given by a similar 
relation to ([3]) replacing f u -i and f u by r£ and Pu ■ This can be done also in order to have minimum values using ( 
i). 



Phyllotaxis with noble number 

Phyllotaxis is characterized by two parameters. One, A, describes the azimuthal (circle) map 6(s + 1) = 6(s) + 2ir\, 
independently of the radial (growth) map r(s). The radial map defines the generative spiral r{6). On a flat substrate, 
the parameter a describes all kinds of parabolic spirals r{9) = a9 a , but also the equiangular spiral r{9) = aexp£> 



(a = co ) and the dense spiral r(9) = aln{l + 8} (a = 0) [17j . On a curved substrate, the spiral is defined by 
projection and an effective growth parameter a > can be defined, with ao = 1/2 for a flat substrate, a > ao and 
a looser spiral for a substrate with positive curvature, and a < ao and a tighter spiral for a substrate with negative 
curvature (our next paper). In this paper, we concentrate on the "ideal" flat case, A = 1/t, ao = 1/2 that yields a 
sequence of concentric circular grain boundaries. 

An infinite phyllotactic structure can be obtained from any < A < 1 that is a noble irrational, A = {qi , q% , . . . ,q u , 1 /r} , 
beginning with an arbitrary rational {qi, (72, ... , q u } and terminating with a golden tail 1/t = {1, 1, . . . , 1, . . .} (contin- 
uous fraction expansion with qi > 1). Now, all rationals between and 1 can be arranged and ordered on a Farey tree, 
and the golden tail oscillates down the tree starting from its defining rational {(71,(72, ■■■ leading to a cascade 



of regular (parastichy) transitions through grain boundaries [3, [la l23l - l26l | . Note that the topological mechanism is 
always the same, articulated through the golden mean 1/t = {1, 1, . . . , 1, . . .}. The actual value of the parameter A 
also oscillates between narrower bounds as the structure becomes larger. This implies also that the measurement of 
A from a given, finite structure is not that precise, and that it is collimated as the structure grows in size. Note that 
there are no exception in nature (ie no false phyllotaxis). Such tree already appears in Koch's thesis and in Rothen 



and Koch II 18[, but the remark that it was isomorphic to Farey is due to Rivier[23j. It was made independently by 
Levitov 24[. Regular and singular transitions were first introduced by Koch and Rothen 12611 . In Levitov, the singular 
transition branches are even disconnected (as an energy flow) from the main Farey tree[27|. 

The corresponding structure has a more or less disorganized central core moulded by the starting rational 
{?!, (72, • . . , qu]-, followed by the well established phyllotactic structure moulded by the golden tail. Obviously, 
the structure must be (much) larger than its core to exhibit it. Note that even for A = 1/t there is an apparent core 
for s < 30 within which successive grain boundaries touch or even overlap. Topologically, the structure is perfect 
down to s = 1. (Table 1). Overlap occurs on strongly negative substrate (our next paper). 

False phyllotaxis exhibit, beyond their core, spider-web structures, generated by A that arc cither rationals, A = 
{<7i, (72, ■ ■ ■ j Qui °o} or Liouville irrationals A = {(71, (72, ■ ■ • , Qu> Q, ■ ■ ■}, with the integer Q 3> qi > 1. On the Farey tree, 
the Liouville irrational goes through a large number < Q of singular transitions. 



CONCLUSIONS 



An infinite number of sites at distance d organized on the nodes of a triangular tiling defined by a hexagonal lattice 
have the translational and rotational symmetries of this lattice. To each nodes it is possible to associate an hexagon 
considered as its Voronoi cell. This tessellation with regular hexagonal cells, obeys the Euler relation F — E + V = 
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where F, E and V are respectively the number of faces, edges and vertices in the tiling. If we consider a limited 
number of points occupying a finite domain limited by a circular boundary, the Euler relation becomes F — E + V = 1 
and non-hexagonal cells appear among the hexagonal cells. The distribution of points in the disk therefore contains 
inherent defects and is no longer invariant under the symmetry operations of classical crystallography. 

The algorithm of phyllotaxis builds a dense organization of points in a disk that optimizes the homogeneity of 
the area associated to each point and the isotropy of the environment of each point in a situation of a circular 
symmetry. To identify the laws specific to this type of organization we have described it using Voronoi and Dclaunay 
decompositions which make apparent the presence of intrinsic defects. These defects are pentagonal and heptagonal 
cells associated two by two in dipoles organized along narrow concentric rings separating larger rings containing 
hexagonal cells continuously deformed. Such an organization can be understood by considering two types of disorder 
in the organization: a metric disorder corresponding to fluctuations of distances between first neighbors, as the 
hexagonal cells are not regular, and a topological disorder corresponding to the presence of pentagonal and heptagonal 
cells associated into dipoles. These two disorders interact to build the organization. Notice that as usual, dislocations 
are called topological defects, but the notion of defect refer to the perfect crystalline structure. In the phyllotaxis 
pentagon-heptagon dipoles are a basic constituent of the structure. 

The role of pentagonal and heptagonal cells dipoles is to introduce rows of new cells needed to maintain the 
homogeneity of area when the radius of the structure increases, as does a dislocation in metallurgy. This organization 
is reminiscent also of grains and grain boundaries of dislocations in metals, but here the grain boundaries are imposed 
by the circular symmetry whereas in metal grain boundaries they result from processing treatments. Finally a very 
specific order is observed in the organization of these dipoles in concentric grain boundaries: the radius of the grain 
boundaries, or the width of grains, follow the Fibonacci series, while the sequence of dipoles along the grain boundaries 
are those of one-dimensional quasicrystalline sequences which are deduced from each other by an inflation-deflation 
rule which leads to a self-similar structure. 

Thus the topological constraint of the circular symmetry introduces an original inflation-deflation symmetry replac- 
ing the translational and rotational symmetries of classical crystallography. Still using the language of metallurgy, 
the Voronoi cell and parastichies could be seen as the repeat unit and the crystallographic lattice planes of this 
unconventional structure. 

Rings of dipoles, the grain boundaries, with their self-similar organization are the fundamental characteristic of 
phyllotaxis. Presently wc are investigating phyllotactic tilings on the sphere and on the hyperbolic plane, always 
keeping a constant area for each site. The geometry of the dipole rings and their succession is the same as on the flat 
disk. All the effect of the Gaussian curvature of the tiled surface only concerns the width of the hexagonal grains. 
So whatever is the underlying geometry grain boundaries are structurally blocked, without any adjustable parameter. 
This is true with a topological but also a metric point of view. 

The organization of the florets of the flowers of many plants is the earliest known example of implementation of the 
algorithm of phyllotaxis, it was indeed at the origin of the development of this domain of studies. As already quoted 
in the introduction recent works on purely physical systems offer new examples [8|, llll |28| . 

All these examples are macroscopic assemblies of elements such as florets, leaves, convection cells, fcrrofluid droplets 
and bubbles. They have in common, not only a circular symmetry, but also the ability of their elements to move with 
respect to each other in order to adjust their distribution in the curse of the growth process. Those symmetry and 
behavior are also presented at the microscopic level by molecular associations of soft condensed matter and biological 
materials. Their molecular interactions allow an internal mobility and their interfacial tension in solution imposes 
the circular symmetry so that they could be relevant of a similar approach. This is why we proposed it for describing 
collagen fibers [li|. If this were so, whereas in the cases of flowers or droplets the growth is generated from the center 
and in that of the convection cells the organization builds itself at a collective transition, this growth of fibrils by 
addition at their surface would provides a third example of the power of the algorithm of phyllotaxis in structural 
studies. 
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